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Abstract 

In this paper we present the results of Monte Carlo simulation studies on the capability of the proposed NEMO 
km 3 telescope to detect high energy neutrinos. We calculated the detector sensitivity to muon neutrinos coming 
from a generic point-like source. We also simulated the lack of atmospheric muons in correspondence to the Moon 
disk in order to determine the detector angular resolution and to check the absolute pointing capability. 
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1. Introduction 

The NEMO Collaboration is conducting an 
R&D activity towards the construction of a 
Mediterranean km 3 neutrino telescope [1]. In this 
work, we present a first study of the expected re- 
sponse of the proposed NEMO km 3 detector to 
neutrinos from point-like sources. In section 2, we 
describe the generation of the atmospheric neu- 
trinos and muons background. We define the de- 
tector sensitivity and the event selection criteria. 
Eventually, we determine the detector sensitivity 
to a muon neutrino flux from a generic point-like 
source. In section 3, we present results on possible 
Moon shadow detection. We then estimate the de- 
tector angular resolution and we check the effects 
of errors in the detector absolute orientation in 
the pointing accuracy. 



Email address: carla.distefano@lns.infn.it (C. 

Distefano). 



The geometry of the km 3 telescope, simulated 
in this work, is the NEMOdh-140 described in 
[2]. The detector response is simulated using the 
ANTARES simulation codes [3], taken into ac- 
count the water optical parameters measured in 
the site of Capo Passero [4] . 



2. Detector sensitivity to point-like sources 

The detector sensitivity is here calculated fol- 
lowing the the Feldman & Cousins approach [5] . In 
particular we define the sensitivity as: 

( &£v\ = 7%)») / ^pA 
\ds„J 90 n s / q ' 

where ~pgo(b) is the 90% c.l. average upper limit 
for an expected background with known mean b 
and (d(p u /de u )o is an arbitrary source spectrum 
predicting a mean signal n s . 
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2.1. Simulation of the atmospheric background 

A sample of 7 • 10 9 atmospheric neutrinos have 
been generated using the ANTARES generation 
code, in the energy range 10 2 10 8 GeV, with a 
spectral index X=2 and a 4tt isotropic angular dis- 
tribution. The events are weighted to to the sum 
of the Bartol flux [6] and of prompt neutrino rqpm 
model [7] flux. Atmospheric muons are generated 
at the detector, applying a weighted generation 
technique. We simulated 5.5 • 10 7 events in the en- 
ergy range 10 2 -=- 10 6 GeV. The events are weighted 
to the Okada parameterization [8] , taking into ac- 
count the depth of the Capo Passero site (D = 
3500 m) and the flux variation inside the detector 
sensitive height. 

2.2. Criteria for atmospheric events rejection 

The used reconstruction algorithm is a robust 
track fitting procedure based on a maximization 
likelihood method. We use, as goodness of fit crite- 
rion, the variable: 
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where log(£)/ Ndof is the log-likelihood per de- 
gree of freedom and N comp is the total number of 
1° compatible solutions [3]. A quality cut A > A cut 
is applied together with other selection criteria: 
the muon must be reconstructed as up-going; the 
event must be reconstructed with a number of hits 
Nfu > Nff*] only events reconstructed in a circu- 
lar sky region centered in the source position and 
having a radius of r^n are considered. The optimal 
values of A cnt , and run are chosen to mini- 

mize the detector sensitivity, taking into account 
both atmospheric neutrino and muon background. 



2.3. Sensitivity to a generic point-like source 

We simulated muons induced by ~ 10 9 neutri- 
nos with energy range 10 2 -j- 10 8 GeV and X=l. 
These events are weighted to the neutrino spec- 
trum (dip v /de v )o = 10~ 7 e~Q eV (GeV -1 cm -2 
s -1 ). As a representative case, the source position 
is chosen at a declination of 5 = —60°. In Table 



1, we report the expected sensitivity for different 
spectral index a, compared to the IceCube tele- 
scope [9]. In Fig. 1 we plot the NEMO sensitivity 
for a = 2 as a function of years of data taking, 
compared to the IceCube sensitivity obtained for 
a 1° search bin [9]. The proposed NEMO detector 
shows a better sensitivity to muon neutrino fluxes, 
reached in a smaller search bin. 



NEMO 


IceCube 


a Acut Nff* r bin #90 (b) e™(dip u /dev)vo 


e%(d(p u /d£ u )90 


1.0 7.6 26 0.4° 2.43 1.9 • 1CT 15 
1.5 7.3 - 0.5° 2.45 2.6 • 1CT 12 
2.0 7.3 - 0.5° 2.45 1.2 • 10~ 9 
2.5 7.3 - 0.6° 2.48 2.3 • KT 7 


2.4- 10- 15 
4.5 • 10- 12 
2.4 • 10 -9 
3.8 • 10- 5 



Table 1 

Sensitivity to a point-like neutrino source at S = —60°, for 
different spectral indexes a and 3 years of data taking, and 
comparison with the IceCube detector [9]. The sensitivity 
spectrum e^(d(p u /de u ) 90 is expressed in GeV a_1 / cm 2 s. 




Years of Data Taking 

Fig. 1. Sensitivity to a neutrino spectrum with a = 2, 
coming from a 5 = —60° declination point-like source and 
comparison with the IceCube detector [9]. 



3. Detection of the Moon shadow 

Since the Moon absorbs cosmic rays, we expect a 
lack of atmospheric muons from the direction of the 
Moon disk (angular radius Rmoou — 0.26°). The 
detection of this muon deficit, commonly called 
Moon shadow, provides a measurement of the de- 
tector angular resolution; besides detecting its po- 
sition in the sky allows us to determine the abso- 
lute orientation of the detector [10]. In this section, 
we present first results of our Monte Carlo simula- 
tions. 
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3.1. Simulation of atmospheric muons 

Atmospheric muons are generated as described 
in section 2.1, we implemented the calculation of 
the Moon position to simulate the muon lack. We 
simulated 1.25 -10 8 events in the energy range 10 2 -j- 
10 6 GeV, restricting the generation in a circular 
window around the Moon position with a radius 
of 10°. The events are then weighted to the Okada 
parameterization. 

3.2. Estimate of the detector angular resolution 



The number of muons reconstructed within 
a distance D from the Moon centre is plotted in 
Fig. 2. The distribution is then fitted using the 
function: 
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obtaining from the fit an angular resolution a = 
0.19° ±0.2°. During the fit procedure, we apply an 
event selection imposing A cut = — 7.6 and NJ 1 ™ = 
20, reaching, in this way, a detection significance 
S'lyr ~ 5.3 in 1 year of data taking. Assuming a re- 
quired significance of 3, the minimum time needed 
to observe the shadow is therefore of the order of 
100 days. 



simulated, introducing a rotation \I/ in the recon- 
structed tracks around the Z axis. In Fig. 3, we 
plot the selected events in a 2D celestial map re- 
ferring to the Moon position, for ^ = 0°, 0.2°, 0.4° 
and 0.6°. For the expected accuracy ^ < 0.2°, 
the shadow is still observable at the Moon posi- 
tion, while considering the pessimistic case \£ > 
0.2°, systematic errors may be corrected. A de- 
tailed study together with the effects of errors in 
the absolute zenithal orientation are still under 
analysis. 




Fig. 3. Celestial map of selected events, referring to the 
Moon position, for ^ = 0° (top-left), 0.2° (top-right), 0.4° 
(bottom-left) and 0.6° (bottom-right). The black circle rep- 
resents the Moon disk. 




D (deg) 

Fig. 2. Selected muon event density versus the angular 
distance from the Moon centre. The line represents the 
result of the fit in Eq.(3), obtained for k = (659±8) deg -2 
and a = 0.19° ±0.2°. 



3.3. Check of the detector pointing accuracy 

Up to now, we assumed to well know the detec- 
tor absolute orientation. The effect of a possible 
error in the absolute azimuthal orientation is here 
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